Effects of the Detection Efficiency on Multiplicity Distributions 



A.H. Tang 1 , G.Wang 2 
1 Brookhaven National Laboratory, Upton, New York 11973 
2 University of California, Los Angeles, California 90095 



m 
o 



X 

=5 
S3 



> 
OS 

m 

O 

m 



x 

S3 



In this paper we investigate how a finite detection efficiency affects three popular multiplicity 
distributions, namely the Poisson, the Binomial and the Negative Binomial distributions. We found 
that a constant detection efficiency does not change the characteristic of a distribution, while a 
variable detection efficiency does. We layout a procedure to study the deviation of moments from 
the baseline distribution due to a variable detection efficiency. 

PACS numbers: 25.75.Ld 



I. INTRODUCTION 

One of the main purposes of relativistic heavy ion col- 
lision experiments is to explore the QCD phase bound- 
ary [l[ , in particular to look for signatures of a first order 
phase transition 0, Q and a critical end point [3, ■ Mo- 
ments of the distributions of conserved quantities, such as 
net-baryons, net-charge and net-strangeness, have been 
argued to be sensitive to the phase transition and the crit- 
ical end point, and are drawing increased attention from 
both experimentalists and theorists [^-[T^j . In the 

study of higher order moments and their derivative quan- 
tities, an abnormal deviation from the baseline distribu- 
tion is usually interpreted as an interesting physics signal. 
In practice, such a deviation is complicated by experi- 
mental effects, such as a finite detection efficiency. In this 
paper, we address how a finite efficiency would change 
three widely used multiplicity distributions, namely, the 
Poisson, the Binomial and the Negative Binomial distri- 
butions. We will discuss the case of a constant efficiency, 
followed by the variable-efficiency case, where we lay- 
out a procedure to investigate how the efficiency affects 
the three multiplicity distributions. The procedure also 
applies to the difference distribution of two multiplicity 
distributions. 



II. MULTIPLICITY DISTRIBUTIONS WITH A 
CONSTANT EFFICIENCY 

A. Poisson Distribution 



g(y) = (i 



(3) 



where e is the probability of seeing a particle, in practice 
less than unit due to the finite acceptance and detection 
efficiency. Without losing generality, below we refer to e 
as the detection efficiency for both sources combined. 
Then Eq. ^ becomes 
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One immediately identifies that the new generating 
function, for an experimental observable with a finite de- 
tection efficiency, still maintains the form of a Poisson 
distribution, with the mean of the distribution reduced 
to Ae. 



B. Binomial Distribution 



f(k;n,p) = (fjp k (l- P y 
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Its probability-generating function is given by 

G(z) = (l-p + pz) n . 
Similarly with a finite detection efficiency, 
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f(k;X) = 



X k e~ x 
k\ ' 



Its probability-generating function is given by 
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We treat observing and not-observing a particle as "de- 
cay" modes of a particle, and apply the cluster decay 
theorem 13[ by replacing z with the generating function 



G(z)=G(g(y)) = [l-p+ p (l- e 
= [1 -pe + pey] n . 
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We have recovered the probability-generating function 
for the Binomial distribution with the replacement of 
p — > p'(= pe). The mean of the new distribution is given 
by p! = fie. The calculation of other quantities under the 
influence of a finite detection efficiency is thus straight- 
forward. For example, 



KG 



cu 
c 2 



1 — 6p + 6p , 
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where k is the kurtosis and Cj are ith order cumulants. 
When taking the detection efficiency into account, one 
simply replaces every p with pe, 



Ca 



kg = — - = 1 — 6pe + 6pV. 

Co 



(9) 



Such knowledge is useful for quantifying the deviation of 
the observable of interest from the original distribution 
due to the finite detection efficiency. 



C. Negative Binomial Distribution 



f(k;r,p)= {^ + r k l \l-pfp r , (10) 

with identities of p = A^ and r = A^-, where /i and a 2 
are the mean and variance, respectively. 

Its probability-generating function has the form of 



G(z) 



1 + z-z 
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(11) 



where p — At — -A- . 

Likewise, in the case of a finite detection efficiency, we 
have 



G(z) = G(g(y)) = 



P 



l-(l-p)(l-e + ey) 
P' 



\-{\-p')y 



(12) 



where p' = —A? — , and r is unchanged. The form 
of the probability-generating function for the Negative 
Binomial distribution is recovered, with p — > p' and 
H — > n (= /ie). Again, other quantities with a finite de- 
tection efficiency can be evaluated with the two simple 
replacements. For example, replacing p with 1 _ p ^_ £ - ) 
everywhere in 
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gives the no 2 for the case with a finite detection efficiency. 



III. MULTIPLICITY DISTRIBUTIONS WITH A 
VARIABLE EFFICIENCY 

Usually the detection efficiency decreases with in- 
creased multiplicity, as the reconstruction of a parti- 
cle becomes more difficult in a lousy environment. In 



this case, the detection efficiency is expressed as a func- 
tion of k, e(k). Now for all the three distributions, the 
probability-generating function can no longer be written 
in a concise form. Instead, we take the general definition 



G(y) = Y,f( k > k 

oo 

= ^/(fc)[l-6(fc) + e(%] fe . (14) 

k=0 

Generally one cannot recover the generating function 
of the same type. That means, a variable efficiency will 
distort the original distribution, unlike the case of a con- 
stant efficiency, where the detector effect will change the 
mean and width of the distribution, but keep the char- 
acteristic shape (as the same type). Nevertheless, with 
e(k) as input, one can still calculate the mean (//) and 
the variance (a 12 ): 



H f = (M)=A 1 , 

a 12 = (M 2 )-(M) 2 

= (M (M - 1)) + (M) - (M) 2 
= A 2 + A 1 -A\, 

where Ai is the factorial moment (M(M — 1) • • 

1)), given by Ai 



(15) 
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For the Poisson distribution 



Ai = 



e(kY 



for the Binomial distribution 
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and for the Negative Binomial distribution 

r (i-p)V e (fc)'. 



= - (fc + r-1)! 
1 ^(fc-*)!(r-l)!' 
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With Eq. ([!?]), ([I8|) and ([15]), A, can be numerically cal- 
culated with known e(fc), and the calculation is no more 
complicated than that for the corresponding distributions 
with the perfect detection. Note that in practice one only 
needs to perform the summation over k to a value that is 
large enough, say, a few a above the mean value, so that 
Ai has little change with further increase of k [14[ . 
The third and fourth central moments are given by 



((M - (M)f) = A x + 2A\ + 2,A 2 - 2>A 1 {A 1 + A 2 ) + 



and 

<(A/-(M}) 4 ) 



A 3 , 
(20) 



A l - 3Aj + 7A 2 + 6Af (Ai + A 2 )+ 

6A 3 - 4Ai(A! + 3A 2 + A 3 ) + A 4 . 

(21) 
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With mean, variance, the third and fourth central mo- 
ments, the first few cumulants can be calculated as usual 

Ci = ((SM)) 

C 2 = ((SM) 2 ) 

C 3 = ((SM) 3 ) 

d = ((SM) 4 )-3((SM) 2 ) 2 , (22) 

where 5M = M — (M). One can further calculate skew- 
ness and kurtosis based on cumulants, which is straight- 
forward and thus is not repeated here. 

Note that although we addressed three specific multi- 
plicity distributions, the procedure discussed in this sec- 
tion can be extended to other multiplicity distributions, 
as long as the factorial moments can be conveniently cal- 
culated. 



IV. DIFFERENCE DISTRIBUTION OF TWO 
MULTIPLICITY DISTRIBUTIONS 

The difference between two independent variables is 
useful for studying the fluctuation of conserved quanti- 
ties, e.g., the net charge and the net baryon number. The 
difference between two variables, each following the Pois- 
sion distribution, is called the Skellam distribution, and 
its probability-generating function is given by: 

G(z; mi,m 2 ) = e -("i+/*)+f i*W*. (23) 

It follows from one of the properties of the probability- 
generating function: for the difference of two independent 
random variables S = X\ — X 2 , the generating function 
is given by Gs(z) — Gx 1 (z)Gx 2 ( z ~ 1 )- The generating 
function for the difference between two Binomial vari- 
ables is 

G(z;n 1 ,p 1 ,n 2 ,p 2 ) = (l-pi+P\z) ni (l-p 2 +p 2 /z) n \ 

(24) 

and the generating function for the difference between 
two Negative Binomial variables is 



which the argument on the detection efficiency still ap- 
plies. This facilitates the calculation of cumulants of the 
three difference-distributions with the finite detection ef- 
ficiency under consideration. For example, for the net 
charge distribution, the additivity of cumulants directly 
gives CAchargc = C + — C_ , where C + and C_ are cumu- 
lants for positively and negatively charged particles, re- 
spectively. The CAchargc with a finite detection efficiency 
can be calculated as long as C+ and C_ are calculated 
independently. 

V. CONCLUSION 

We have shown that for the Poisson, Binomial and 
Negative Binomial distributions, a constant efficiency 
will modify the mean and the width of the original dis- 
tribution, but it does not change the distribution type. 
With a known constant efficiency, the original distribu- 
tion can be completely reconstructed from the measured 
one, and vice versa. However, a variable efficiency will 
distort the original distribution. In this case the origi- 
nal distribution is difficult to be recovered. Nevertheless, 
one can still study how a finite, variable detection ef- 
ficiency changes the original distribution. With known 
form of e(k), the deviation of moments and their deriva- 
tive quantities from the baseline distributions, can be es- 
timated following the procedure presented in this paper. 
The procedure applies also to the difference distribution 
of two independent distributions. Knowledge obtained 
in this work is useful for the investigation of how the 
observable of interest deviates from the baseline distri- 
bution when a finite detection efficiency is present. This 
will help avoid the mis-interpretation of certain observ- 
ables as signals of the possible phase transition and/or 
the critical end point. 



Acknowledgments 



q, ) — ( Pi \ ( P^ ] The work for the case of a variable detection efficiency 

' ' ' ' \1 — (1 — pi)z J \1 — (1 — p 2 )/z J was stimulated by a discussion with J. Dunlop and P. 

(25) Sorensen. A. Tang was supported by the US Department 

When we take into account the finite detection efficiency, of Energy under Grants No. DE-AC02-98CHI0886 and 

none of the three generating functions above can re- No. DE-FG02-89ER4053I . G. Wang was supported by 

cover the form of the same type. Fortunately, they de- the US Department of Energy under Grants No. DE- 

scribe the difference between two quantities, to both of FG02-88ER40424. 



[1] I. Arsene et al. Nucl. Phys. A 757, 1 (2005); 
B. B. Back et al. Nucl. Phys. A 757, 28 (2005) 
J. Adams et al. Nucl. Phys. A 757, 102 (2005); 
K. Adcox et al. Nucl. Phys. A 757, 184 (2005); 
S. Chatrchyan et al. JHEP 1205, 063 (2012); 



B. Abelev et al. Phys. Rev. Lett. 109, 072301 (2012). 
[2] S. Ejiri, Phys. Rev. D 78, 074507 (2008). 
[3] E. S. Bowman and J. I. Kapusta, Phys. Rev C 79, 015202 

(2009). 

[4] M. A. Stephanov, Prog. Theor. Phys. Suppl. 153, 139 



4 



(2004); 

[5] Z. Fodor and S. D. Katz, JHEP 0404, 050 (2004). 
[6] M. M. Aggarwal et al. Phys. Rev. Lett. 105, 022302 
(2010). 

[7] D. McDonald [for the STAR Collaboration], 

larXiv:1210.7023l (2012) 
[8] J.T. Mitchell [for the PHENIX Collaboration], 

I arXiv:1211. 61391 (2012) 
[9] M. A. Stephanov, Phys. Rev. Lett. 102, 032301 (2009). 
[10] M. Asakawa, S. Ejiri and M. Kitazawa, Phys. Rev. Lett. 



103, 262301 (2009). 

[11] C. Athanasiou et al, Phys. Rev. D 82, 074008 (2010). 

[12] M. A. Stephanov, Phys. Rev. Lett. 107, 052301 (2011). 

[13] J. Pumplin, Phys. Rev. D 50, 6811 (1994). 

[14] To avoid the problem of the overflowing with large nu- 
merical numbers, one may also calculate, for example, 

oo , OO , 

e- X E T^eikY as exp[log(£ T^AW) - A]. 

k — i k — i 



